PROCEEDINGS 

OF SCIENCE 



Higher loop renormalization of fermion bilinear 
operators 



A. Skouroupathis and H. Panagopoulos 

University of Cyprus, Department of Physics, Cyprus 



E-mail: php4as01@ucy.ac.cy haris@ucy.ac.cy 



We compute the two-loop renormalization functions, in the RI ' scheme, of local bilinear quark 
operators \j/T\j/, where F denotes the Scalar and Pseudoscalar Dirac matrices, in the lattice for- 
mulation of QCD. We consider both the flavor non-singlet and singlet operators; the latter, in the 
scalar case, leads directly to the two-loop fermion mass renormalization, Z m . 
As a prerequisite for the above, we also compute the quark field renormalization, Zy, up to two 
loops. 

We use the clover action for fermions and the Wilson action for gluons. Our results are given 
as a polynomial in c$\y, in terms of both the renormalized and bare coupling constant, in the 
renormalized Feynman gauge. We also confirm the 1-loop renormalization functions, for generic 
gauge. 
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1. Introduction 

Studies of hadronic properties using the lattice formulation of QCD rely on the computation 
of matrix elements and correlation functions of composite operators, made out of quark fields. 
A whole variety of such operators has been considered and studied in numerical simulations, in- 
cluding local and extended bilinears, and four-fermi operators. A proper renormalization of these 
operators is most often indispensable for the extraction of physical results from the lattice. 

In this work we study the renormalization of fermion bilinears 6 = yTy on the lattice, where 
r = 1 , 75 . We consider both flavor singlet and nonsinglet operators. The cases r = y^ , 75 7^ , 75 cfy v , 
will be presented in a sequel to this work. In order to obtain the renormalization functions of 
fermion bilinears we also compute the quark field renormalization, Zy, as a prerequisite. 

We employ the standard Wilson action for gluons and clover-improved Wilson fermions. The 
number of quark flavors Nf, the number of colors N c and the clover coefficient csw are kept as free 
parameters. Our two-loop calculations have been performed in the bare and in the renormalized 
Feynman gauge. For 1-loop quantities, the gauge parameter is allowed to take arbitrary values. 

The main results presented in this work are the following 2-loop bare Green's functions (am- 
putated, one-particle irreducible (1PI)): • Fermion self-energy: Ly(q,a L ), • 2-pt function of the 
scalar \ff\ff : Y^{qa L ), • 2-pt function of the pseudoscalar xjfysyf : T.p(qa L ) (a L : lattice spacing, 
q : external momentum) 

In general, one can use bare Green's functions to construct Zj , the renormalization function 
for operator G, computed within a regularization X and renormalized in a scheme Y. We employ 

L RI 1 L RI 1 L RI 1 

the RI scheme to compute the various renormalization functions: Zy , Z s ' , Z p ' . The corre- 
sponding quantities in the MS scheme: Zy MS , Z$ ,MS , Zp' MS , are presented in the longer write-up 
[§■ 

The flavor singlet scalar renormalization function is equal to the fermion mass multiplicative 
renormalization, Z m , which is an essential ingredient in computing quark masses. 

Finally, as one of several checks on our results, we construct the 2-loop renormalized Green's 
functions in RI': L^'(q,p.) (6 = y/,S, P), as well as their counterparts in ~MS: L^ s (q,p,). The 
values of all these functions, computed on the lattice, coincide with values computed in dimensional 
regularization (as can be inferred, e.g., from [0]). 

The present work is the first two-loop computation of the renormalization of fermion bilinears 
on the lattice. There have been made several attempts to estimate Z@ non-perturbatively; see Ref.[jl|] 
for a list of recent references. Some results have also been obtained using stochastic perturbation 
theory [[§]. A related computation, regarding the fermion mass renormalization Z m with staggered 
fermions can be found in [Q]. 

2. Formulation of the problem 

Lattice action: We will make use of the Wilson formulation of the QCD action on the lattice, with 
the addition of the clover (SW) term for fermions. The clover coefficient csw in the Langrangian 
is treated here as a free parameter; The fermion mass m a in the Langrangian is a free parameter. 
However, since we will be using mass independent renormalization schemes, all renormalization 
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functions which we will be calculating, must be evaluated at vanishing renormalized mass, that is, 
when m is set equal to the critical value m cr : m Q — ► m cr = + &{g%). 

Definition of renormalized quantities: As a prerequisite, we will need the renormalization func- 
tions for the gluon, ghost and fermion fields (A" , c a , y), and for the coupling constant g and gauge 
parameter a, defined as follows: 

A£ = V^A£ , c a Q = VZc a , Yo = v %¥ , go = H e Z g g , a = Z- l Z A a (2.1) 

The scale /I enters the relation between g a and g only in dimensional regularization (D = 4 — 2e 
dimensions). We will need Za,Z c , Z a andZg to 1 loop and Z v to 2 loops. 

Definition of the RI' scheme: This renormalization scheme [{7|, |8[ ^] is more immediate for a lattice 
regularized theory. It is defined by imposing a set of normalization conditions on matrix elements 
at a scale p., where (just as in the MS scheme): jl = \i (An / e^) 1 1 2 {Je is the Euler constant) 
In Euclidean space, the fermion self energy L^(q,a L ) = i4 + m Q + &{g 2 ) is renormalized by: 



lim 

a, -^0 



Z£*> L /X)tr (L L ¥ (q,a L )4)/(4iq 2 



1 (2.2) 



The trace here is over Dirac indices; a Kronecker delta in color and in flavor indices has been 
factored out of the definition of E£. Similar conditions hold for Z c (extracted from the ghost 
self energy) and Za, Z a (extracted from the gluon propagator). We have checked explicitly that 

L RI 1 

Za = 1 up to one loop, in agreement with the continuum. 

For consistency with the Slavnov-Taylor identities, Z g in the RI' scheme is defined as in the 
MS scheme. In dimensional regularization (DR) this is achieved by tuning the value of Z g in such 
a way as to absorb only the poles in e which appear in the gluon-fermion-antifermion 1PI vertex 
function GAyy (together with matching powers of ln(47r) — }fe); this leads to a result for 
which is finite but not unity. Alternatively, a similar procedure can be performed on the gluon- 
ghost-antighost vertex Gasc, leading to exactly the same value for Z g . 

The corresponding renormalization condition can be applied on the lattice vertex, G\^, re- 
quiring that the renormalized expression G^^ is the same as that stemming from DR, and similarly 

J L RI* j j 

for G AEc . We have calculated Z g , using either one of G ASc or G A ^, and have verified that the 
same result is obtained. 

Renormalization of fermion bilinears: The renormalization of lattice operators &y = i s 
defined by: Gf' = z[ RI '(a L jEx) & To 

L RI 1 

For the scalar (S) and pseudoscalar (P) operators, Z r ' can be obtained through the corre- 
sponding bare 2-point functions ^(qa^) (amputated, 1PI) on the lattice, in the following way: 



lim 

a, ^0 L 



ryL,Rl' ryL.Rl' V L/ \ 



1 , lim 

a, ^0 L 



ZL,Rl' ryL,Rl' j 



7s (2-3) 

L.RI' 



Once Z r ' have been calculated, one may convert them to the MS scheme (see Ref.[[T]]). 



3. Computation and Results 

There are 28 Feynman diagrams contributing to the fermion self-energy L^(q,a L ) at 1 and 2 
loops, and 21 diagrams contributing to Lg(qa L ), Lp(qa L ). For flavor singlet bilinears, there are 4 
extra diagrams, shown in Fig. 1 , in which the operator insertion occurs inside a closed fermion loop. 
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The evaluation and algebraic manipulation of Feynman diagrams, leading to a code for nu- 
merical loop integration, is performed automatically using our software for Lattice Perturbation 
Theory, written in Mathematica. The most laborious aspect of the procedure is the extraction of the 
dependence on the external momentum q. This is a delicate task at two loops; for this purpose, we 
cast algebraic expressions (typically involving thousands of summands) into terms which can be 
naively Taylor expanded in q to the required order, plus a smaller set of terms containing superficial 
divergences and/or subdivergences. The latter can be evaluated by analytical continuation to D > 4 
dimensions, and splitting each expression into a UV-fmite part (which can thus be calculated in the 
continuum), and a part which is polynomial in q. 

Some of the diagrams contributing to Ly(q,a L ), L^(qa L ) and Lp(qa L ) are infrared divergent 
when considered separately, and thus must be grouped together in order to give finite results, for ex- 
ample, diagrams (1-2), (3-4) in Fig. 1. As mentioned before, all calculations should be performed at 
vanishing renormalized mass; this can be achieved by working with massless fermion propagators, 
provided an appropriate fermion mass counterterm is introduced on 1-loop diagrams. 







12 3 4 

FIG. 1. Extra two-loop diagrams contributing to Zg singlet A cross denotes an insertion of a flavor singlet 
operator. Wavy (solid) lines represent gluons (fermions) 

All two-loop diagrams have been calculated in the bare Feynman gauge (a D = 1). One-loop 
diagrams have been calculated for generic values of a ; this allows us to convert our two-loop 
results to the renormalized Feynman gauge (a R // = 1 or = 1). 

Numerical loop integration was carried out by our "integrator" program, a metacode written 
in Mathematica, for converting lengthy integrands into efficient Fortran code. Two-loop integrals 
were evaluated for lattices of size up to L = 40; the results were then extrapolated to L — > oo. The 
extrapolation error can be estimated quite accurately (see, e.g. Ref. |[T^]), given that L-dependence 
of results can only span a restricted set of functional forms. 

Our 1-loop results for Zy, Z$ and Zp confirm existing results ^ (there is, however, a difference 



in Zp' MS , due to an extra finite renormalization factor, Z5, which is required in MS [|11|]). 
Two-loop results: The evaluation of all Feynman diagrams leads directly to the bare Green's 
functions E^, £| and Lp. These, in turn, can be converted to the corresponding renormalization 



L Y L Y L Y 

functions Zy , Z s and Zp (Y = RI or MS), via Eqs.([2.2[) and (p.3|). To this end, we need the 



L RI ! 

one-loop expression for Z A .To express our results in terms of the renormalized g, we also need 

L RI 1 L RI 1 L RI 1 

the one-loop expression for Z g ' . Our 1-loop computation of Z A ' and Z g ' is in agreement 
with older references (see e.g. [g]). We present below Zy , Z s ' and Zp to two loops in 
the renormalized Feynman gauge a s/ / = 1. Bare Green's functions can be easily recovered from 
the corresponding Z's. The errors result from the L — > 00 extrapolation. We employ a standard 
normalization for the algebra generators: ti(T a T b ) = 8 ab /2. (a = g^/\67i 2 , c F = (Af? — 1)/ (2N C )) 
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yL.Rl i . 

1 -+- Cl Cf 



+a 2 c F 



ln(^jU 2 ) + 11.852404288(5) - 2.248868528(3) c sw - 1 .397267102(5) c| w 

N, 



2 3 3 3 

+ ln(a 2 i u 2 )(- 6.363 17446(8) + 0.79694523 (2) N/ c sw - 4.712691443(4) AVc sw 
+49.83082185(5) c F - 2.24886861(7)c F c sw - 1.39726705(1) c F c| w + 29.03029398(4)^ 
+N f ( - 7.838(2) + 1 . 153(1) c sw + 3.202(3) c| w + 6.2477(6) c^ w + 4.0232(6) c 4 sw 



ZL.Rl' i i 
o = 1 + a Q c F 



+a 2 c F 



+c F (^505.39(1) -58.210(9) c sw +20.405(5) cl w + 18.8431(8) 4 W + 4.2793(2)c£ w 
+N C ( -20.59(1) -3. 190(5) c sw -23. 107(6) c 2 sw - 5.7234(5) c| w - 0.7938(1) c£ w ) 

3 ln^ 2 /* 2 ) - 17.9524103(1) -7.7379159(3) c sw + 1.38038065(4) c| w 

+ ln(a 2 jU 2 ) ( - 8.1721694(5)^/ + 2.3908354(3) A^/c sw - 14. 13807433 (4) N f c 2 sw 
+66.0780218(9) c F - 23.213749(2) c F c sw + 4.141 1425(3)c F c sw + 55. 7975008 (9) N, 
+N f (24.003(3) + 11.878(5) c sw +25.59(1) cl w + 22.078(3) clw - 6.1807(8) c£ w 
+c F ( - 602.35(6) + 66.80(7) c sw + 75.42(5) c sw - 27.759(4) c 3 sw - 2.688( 1) c$ w 

+N C ( - 38.16(4) - 120.26(5) c sw - 16.18(3) c| w + 12.576(3) c| w + 1.0175(8) c| w 



(3.1) 



(3.2) 



Z p = 1 + a Cf 



3 ln(a ,2 /i 2 ) — 27.5954414(1) + 2.248868528(3) csw — 2.03601561(4) c$w 



+a Q c F 



In 2 (a 2 h fi 2 )(^c F +N f -^N^ 



+ ln(a 2 p 2 )(^-8.1721694(4)Ar/ + 2.39083540(6)A^/csw- 14.13807433(4)^V/c| w 
+37.1489292(7)c F +6.746606(l)c F c sw - 6. 1080465(3) c F c 2 sw + 55. 7975008 (7 )N, 
+N f (38.231 (3) - 7.672(5) c sw + 55.32(1) c 2 w - 7.049(3) c 3 sw +4.7469(8) c| w ) 
+c F ( - 876.98(4) + 85.80(2) c S w + 37.37(4) c sw + 19.974(3) c| w + 2.873( 1) c 4 sw 

+N C ( - 104.35(3) -38.70(2) c sw - 13.93(3) c| w -4.429(2) c 3 sw - 1.2898(7) c| w 



(3.3) 



All expressions reported thus far for Z$ and Zp refer to flavor non singlet operators. In the case 
of Zp, all diagrams of Fig. 1 vanish, so that singlet and non singlet results coincide. Z$ on the other 
hand, receives an additional finite contribution (which is the same for the MS scheme as well): 



7 L,RI 
^S, singlet 



Z L S R1 ' + a 2 c F N f { - 107.76(1) + 82.27(2)c sw - 29.727(4) c 2 sw 
+ 3.4400(7) cl w + 2.2758(4) c 4 sw ) 



(3-4) 
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In Figs. 2, 3 and 4 we plot Zy RI , Z^' RI and 
Zp , respectively, as a function of csw- For 
definiteness, we have set N c = 3, p. = l/a L 
and /3 G = 2N c /gl = 6.0. Our results up to two 
loops for each Z are shown for both Nf = 
and Nf = 2, and compared to the correspond- 
ing one-loop results. Furthermore, in the 
scalar case, we also present the two-loop re- 
sult for the flavor singlet operator. 



1.14 

1.12 

1.10 

, 1.08 
5 

J > 
N 1.06 

1.04 



1.02 



1.00 
0.0 




FIG. 2. Z^ RI '{a L fi) vs. c sw (N c = 3, p = l/a L , 
j3 = 6.0). Results up to 2 loops are shown for 
Nf = (dashed line) and Nf = 2 (solid line); 
one-loop results are plotted with a dotted line. 
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FIGS. 3, 4. Zg' m ' (a L fi), Zp Kl (a L p.) vs. c sw (N- = 3, p, = l/a L , j3 Q = 6.0). Results up to 2 loops, for the 
flavor non-singlet, are shown for Nf = (dashed line) and Nf = 2 (solid line); 2-loop results for the scalar 
flavor singlet, for Nf = 2, are plotted with a dash-dotted line; one-loop results are plotted with a dotted line. 
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4. Discussion 



As can be seen from Figs.2-4, all 2-loop renormalization functions differ from 1-loop values 
in a significant way; this difference should be taken into account in MC simulations, in order to 
reduce systematic error. At the same time, 2-loop contributions are consistently smaller than 1-loop 
contributions, indicating that the (asymptotic) perturbative series are under control. 

The dependence on the clover parameter c$w is also quite pronounced. In the present work, 
csw was left as a free parameter; its optimal value, as dictated by &(a L ) improvement, has been 
estimated both non-perturbatively and perturbatively (to 1-loop) in the literature. 

Our results regard both the flavor nonsinglet and singlet operators. For the pseudoscalar op- 
erator, these cases coincide, just as in dimensional regularization. The scalar, on the other hand, 
receives an additional finite (a L jl independent) contribution in the flavor singlet case. Zs^ingUt is 
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seen to be equal to the fermion mass renormalization Z m , which is essential in the computation of 
quark masses. We note also that the scalar and pseudoscalar factors, necessary to convert flavor 
singlets from the RI' to the MS scheme, stay the same as their non-singlet counterparts. 

The 2-loop computation of the renormalization functions for the Vector, Axial and Tensor 
bilinears is work currently in progress. 

Besides the strictly local definitions of fermion bilinears, xj/Tyf, one can consider a family 
of more extended operators, with the same classical continuum limit, as dictated by &(a L ) im- 
provement. The renormalization of these extended operators involves more Feynman diagrams; 
however, the computation is actually less cumbersome, since all additional contributions are now 
free of superficial divergences. We will report the results of this computation in a future work. 

Acknowledgements: We would like to thank J. A. Gracey and S. A. Larin for private communi- 
cation regarding their continuum results. 
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